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We study degeneracies between parameters in some of the widely used parametrized modified 
gravity models. We investigate how different observables from a future photometric weak lensing 
survey such as LSST, correlate the effects of these parameters and to what extent the degeneracies 
are broken. We also study the impact of other degenerate effects, namely massive neutrinos and 
some of the weak lensing systematics, on the correlations. 



I. INTRODUCTION 

The idea of modified gravity (MG) has been pursued as 
an alternative to dark energy (DE) in explaining the ob- 
served cosmic acceleration. Most of the proposed DE/MG 
theories predict more or less the same background evolu- 
tion for the universe so that distance measurement data, 
such as CMB and Supernovae (SNe), can not effectively 
discriminate between them. On the other hand, these the- 
ories have different predictions for the evolution of struc- 
ture [1-5] that can be measured to a high precision by 
the ongoing and upcoming large-scale structure surveys 
[6-11]. 

The growth of structure can be tested in a model- 
independent way by parametrizing the evolution equa- 
tions of the cosmological perturbations. At the linear 
level, it is achieved by introducing two scale- and time- 
dependent functions (MG functions) in the equations re- 
lating the scalar metric perturbations, $ and VP, to each 
other and to the matter perturbations [6, 7, 12, 13, 15- 
20]. 

One of the goals of taking such a phenomenological ap- 
proach is to detect possible deviations of perturbations 
from their default evolution, predicted by ACDM +GR, 
in an efficient and model-independent way. Therefore, it 
is important to study the degeneracies between the pa- 
rameters in these models and their impact on the pa- 
rameter constraints. One can investigate the scales and 
rcdshifts where these degeneracies are minimum, how dif- 
ferent observables can break these degeneracies, or what 
combinations of data do it best. This can help to come up 
with an optimum set of parameters for extracting infor- 
mation from data and constraining the MG parameters. 

In this paper we work with two sets of MG functions. 
The first set parametrizes the modified Poisson equation, 
fj,(k, a), and the ratio between the two Newtonian poten- 
tials, j(k, a) [6, 12, 17]. In the second case, we use a) 
and also parametrize the modifications to the weak lens- 
ing (WL) potential by S(fc, a) [12-14]. We consider the 
approximate parametric form of these functions for some 
of the widely used MG theories in the literature and 
quantify how different data sets can improve on degener- 
acy breaking of the parameters. We use CMB tempera- 
ture and polarization data from the Planck satellite [21], 
SNe from an Euclid-like survey [22] and large scale struc- 



ture data from a photometric WL survey such as Large 
Synoptic Survey Telescope (LSST) [23] . We also consider 
pixelations of (fi, 7) as a more model-independent treat- 
ment of these functions. With the pixels of fi and 7, one 
is able to study the degeneracies at particular redshifts 
and scales of interest. 

There arc different physical processes that lead to de- 
generate effects with modified gravity on the cosmological 
observables. For example, massive neutrinos can mod- 
ify the growth of structure on their free-streaming scale. 
This modification can, in principle, be degenerate with 
the effects of MG on the overlapping scales and redshifts 
[24-26]. We also consider the impact of some of the sys- 
tematic effects expected in WL surveys. We use [27, 28] 
where three sources of systematics are modeled for fu- 
ture tomographic surveys: photo-z errors, additive and 
multiplicative errors due to the uncertainty of the point 
spread function (PSF) measurements. 

In Section II, we describe the formalism, MG mod- 
els, observables and experiments we use. In Section III, 
we study the degeneracies between the MG and other pa- 
rameters of our chosen models where different parametric 
functional forms of the MG functions are used. From a 
combination of the upcoming CMB, SNe, galaxy number 
counts (GC) and WL data, we show that parametriza- 
tions with small number of MG parameters can not cap- 
ture a wide range of possible MG features in the data 
while too many parameters mainly introduces degenera- 
cies with no extra benefit. In Section IV, we investigate 
how the additional effects (mass of neutrinos and WL sys- 
tematics) will change the results. The inclusion of such 
degenerate effects would require other complementary 
data to reduce the correlations between the parameters. 
We summarize in Section V. 



II. THE FORMALISM 

A. Parametrization of the linear perturbations 

The line element for a linearly perturbed FRW universe 
can be written in the Newtonian gauge as 

ds 2 = -(1 + 2*) dr 2 + a 2 {t) (1 - 2$) dx 2 , (1) 
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where '5 and $ are the scalar perturbations of the metric, 
r the conformal time and x the spatial coordinate. 

At the level of linear perturbations, in addition to 
and one needs to evolve the matter density contrast 
5 = bpj p and the velocity field perturbation v. Two equa- 
tions are obtained from the conservation of energy mo- 
mentum tensor T»» = [29] 



-(1 + w) 




(2) 



k 2 8-k z a, (3) 



where 9 = kv, w is the equation of state parameter for 
each fluid component and er is the anisotropic stress term 



defined as (p + P)a 



- [kikj 



i<5 



with TTij being 



the traceless component of the energy-momentum tensor. 

The other two equations required to close the system 
are provided by the gravity theory of the model. In the 
case of general relativity (GR), the equations are 

fc 2 $ = -47rGa 2 ^pA, (4) 
fc 2 ($-*) = 12irGa 2 ^2(p + P)<r, (5) 

where pA = p5 + 3(aH/k)v is the comoving density per- 
turbation and the sum is over all the species present in 
the universe. 

In the parametrized framework, Eqs. (4) and (5), are 
modified by introducing the MG functions to accommo- 
date possible deviations from GR (+ACDM) 

fc 2 * = -fx(k, a)AirGa 2 ^{pA + 3(p + P)a} , (6) 
fc 2 [$ - 7(fc, a)*] = fj,(k, a)127rGa 2 + P)a , (7) 

where at late times the shear term, a, is negligible and 
cold dark matter (CDM) perturbations are dominant so 
that the above equations are simplified. The reason for 
using the scalar potential '5 in Eq. (6), as opposed to 
$, is that the growth of matter inhomogeneities is more 
directly related to ^. To see this note that the growth of 
matter on sub-horizon scales can be derived by combining 
Eqs. (2) and (3) in the limit of k/H -» oo 



s + ns + k 2 ^ = 0, 



(8) 



where H = - and time derivative is w.r.t. conformal time. 

a 

In this approach, p(a, k) and ~/(a, k) are generic func- 
tions that parametrize solutions in alternative gravity 
theories and depend on the choice of the initial con- 
ditions. In some theories, one can calculate their ap- 
proximate form analytically in the quasi-static approx- 
imation. For instance, perturbation equations have been 
derived for some classes of scalar-tensor theories in 
[2, 3, 7, 17, 30-32], and for the DGP theory and its higher- 
dimensional extensions in [1, 33-37]. 



B. Experiments and observables 

The data considered here include CMB temperature 
and polarization (T and E) , SNe observations, weak lens- 
ing (WL) shear of distant galaxies, galaxy number counts 
(GC), and their cross-correlations. 



1. CMB and SNe la data 

We assume the expected CMB data from the Planck 
mission [21] of the European Space Agency (ESA) using 
the same parameters as in [7, 38]. In addition, to bet- 
ter constrain the background expansion parameters, we 
include simulated SNe luminosity data for a Euclid-like 
survey [22]. We generate 4012 data points randomly dis- 
tributed in 14 redshift bins from z = 0.15 to z = 1.55. 
The nuisance parameter M 1 is treated as an undeter- 
mined parameter in the analysis. 



2. Large scale structure data 

We assume the GC and WL data by the Large Synoptic 
Survey Telescope (LSST) [23]. LSST is a proposed large 
aperture, ground-based, wide field survey telescope. It is 
expected to cover up to half of the sky and catalogue 
several billion galaxies out to redshift z ~ 3. For LSST, 
we adopt parameters from the recent review paper by the 
LSST [401. 



3. Observables 

We consider observables from the above surveys and 
use the angular power spectra from auto-correlation and 
cross-correlation of these observables for our Fisher anal- 
ysis. The power spectra are calculated using MGCAMB 
code [7, 41] where Eqs. (6) and (7) are implemented. 
With LSST, we consider 10 and 6 tomographic redshift 
bins for GC and WL observations, respectively. 

The details of the implementations are described in 
[381. 



C. MG models 

We choose to work with four models which include 
one, three, five and twelve MG parameters, respectively. 
Table I shows the model parameters that we use in our 
Fisher analysis and the order in which they are sorted. 
In all of the models, we assume spatially flat geometry 
and include the main cosmological parameters (#1-5 and 
A s ). When GC data is included, we assume that the bias 
is scale-independent on the linear scales and introduce 
10 constant bias parameters, one for each photometric 
redshift bin surveyed by LSST. w{z) is either binned or 
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TABLE I. List of the parameters included in our models (I - IV), with the numbers showing the order of their appearance in 
each model. 



considered as a constant and is included in the Fisher 
analysis . 



1. 1-parameter MG model: f(R) theory 

f(R) models can be tuned to reproduce any back- 
ground expansion history, and the remaining relevant 
quantity is the squared Compton wavelength of the new 
scalar degree of freedom fn = df /dR mediating the fifth 
force. In units of the Hubble length squared it is given 
by [2, 42] 



B 



f RR dR (dhxH 



1 + fn d In a \ d In a 



(9) 



where R is the background Ricci scalar. Thus, for a fixed 
background expansion history, different f(R) models can 
be parametrized by the parameter Bq, which is the value 
of B today. It was suggested in [9] that for B « 1 the 
large scale growth in f(R) models can be modeled as: 



fj,(k,a) 



1 



\BnC 



1 + B C 



7(a, k) 



1 



\BnC 



1 



^B C ' 



(10) 



where C = (cka 2 ) 2 /2Hq. This one-parameter model gives 
a good approximation of f(R) theories in the quasi-static 
limit [3, 9]. We choose a fiducial value of Bq ~ 10~ 6 in 
accordance with the forecasted error on this parameter 
from LSST [381. 



where the lower equation governs the perturbations of the 
relativistic particles through light deflection and ISW ef- 
fect, and a) accounts for deviations from the default 
evolution in the ACDM+GR model. It has been shown 
that the following choice can cover a wide range of scalar- 
tensor theories available [43]: 



f i = l 



ca s {k/H ) n 
l + 3ca s {k/H ) n ' 



1 . 



(13) 



It is easy to show that Eqs. (13) can be described in terms 
of (fJ.,j) as 



2 
/<■ 



1 = 1 



ca s (k/H ) n 
l + 3ca s (k/H ) r 

6ca s {k/H Q ) n 
1 + ca s (k/H ) n 



(14) 



We work with the fiducial values of s — 4 and n = 2 
and c = 0.002. The value for c is taken from the 95% 
C.L. constraints on the parameter from the current large 
scale structure data [43]. 



3. 5-parameter MG model: a general BZ parametrization 

For a wide class of scalar-tensor theories, (/i, 7) can be 
described by the following approximate expressions (BZ 
parametrization) [6, 17] 



2. 3-parameter MG model: scalar-tensor theory 

An alternative parametrization to Eqs. (6,7) can be 
written as [8, 13] 

= -AirGa 2 fi(k, a)pA, (11) 
fc 2 ($ + *) = -87rGa 2 £(fc, a) P A, (12) 



1 We include the SNe la nuisance parameter in all these models 
and then marginalize over it. 



. 1 + ftAffcV 

" (a '* ) = 1 + X 2 k 2 a* 



7(0, fc) = 



1 



- p 2 \ 2 k 2 a' 
1 + X 2 k 2 a s 



(15) 



where /3 2l Ai, A2 and s can have different values in 
different MG theories. The BZ parametrization provides 
more flexibility in capturing MG features. We choose a 
set of fiducial parameter values corresponding to a f(R)- 
class model [7]: ft = 4/3, (3 2 = 1/2, s = 4, X 2 = ^X 2 = 
1000 Mpc 2 . 
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4- 12-parameter MG model: model-independent pixelization 

We also consider the less model-dependent case where 
fi and 7 are "pixelized" in z and k space. We choose to 
have three bins in the linear k range (10 -5 — 0.2 h/Mpc) 
and for each k bin, we have a high-z bin (3 < z < 30) 
and a low-z bin (0 < z < 3). In total, there are twelve 
MG pixels, six for p, (# 8-13) and six for 7 (# 14-19). 
w(z) is also binned to one low-z bin (#6) and one high-z 
one (#7). 



III. DEGENERACY OF PARAMETERS 

In this section, we study the correlations between the 
MG and other parameters of the models described in II C. 
For each MG parameter, we show two correlation plots: 

• Un-normalized correlation CV, , where C is the co- 
variance matrix of all parameters, computed by in- 
verting the corresponding Fisher matrix. 

• normalized correlation: C%j = C,j / 'y '{Cn x Cjj). 

The reason for showing the un-normalized correlations 
is that it is hard to judge about the importance of cor- 
relation between two parameters by just looking at the 
normalized correlation values. For example, a parame- 
ter might be highly correlated with the others while be- 
ing very well-constrained, which means that it is well- 
measured anyway. For the un-normalized case, we plot 
the absolute values of correlations on a semi- log plot. 
Also we show the U Cu x Cjj values where i runs over 
all the parameters and j is the index of the MG parame- 
ter we are plotting the correlations for. These values will 
be used as a reference and help us estimate the relative 
degree of correlations between parameters. 

We start from the combination of the CMB and SNe 
la data. We then add the GC and WL data in two 
steps and study how including them breaks the de- 
generacies. Finally, we add all the cross-correlations of 
CMBxCG, CMBxWL and WLxCG. As mentioned, 
CMB+SNe data are always included and mainly con- 
strain the background parameters. For simplicity, we do 
not show CMB+SNe and CMB+SNe+WL combinations 
in the un-normalized plot. 

With CMB+SNe data only, one measures the cosmo- 
logical distances and thus, the background parameters. 
To describe this, we note that the comoving distance to 
a redshift z is defined as 

f z dz 

Jo h q [n m (i + zf + n r (i + z) 4 + n DE (z)] 1/2 ' 

(16) 

where £l m , &r and £Ide(z) are the matter, radiation and 
dark energy fractional densities. Varying the background 
parameters changes the comoving distances to observed 
phenomena such as recombination. This interferes with 



the growth of perturbations by changing the time avail- 
able for perturbations to grow, which can be degenerate 
with changing the growth rate. 

The growth of the matter inhomogeneities is probed 
by large scale structure surveys through estimation of 
the scalar potential in Eq. (8). With the GC data, one 
can measure the density fluctuations of galaxies (p g 5 g ) 
which is related to the total matter density fluctuations, 
Pm&m-, via the bias parameter 

S CJ = bS m . (17) 

In other words, the potential '5 is estimated from the 
matter density fluctuations via Eq. (6) up to a (linear) 
bias factor. Since baryons fall into the same gravitational 
potential as CDM, we can write 

5 g +H5 g <x^6 g . (18) 

Having 5 g measured by GC observations, the growth of 
matter perturbations is directly related to p, bias and 
the Hubble parameter. This demonstrates the source of 
large degeneracies between the MG parameters, bias and 
the background parameters. 

WL data probes the <I> + 5' combination. With a similar 
argument as above, one can see that the evolution of the 
WL potential is determined by p and £ (Eqs. (11) and 
(12)) and hence, a combination of p and 7. Here, the 
bias parameter is not involved and the growth of matter 
is directly dependent on p m 5 m . With WL data alone, 
there is a high degeneracy between the MG parameters 
and n b h 2 , tt c h 2 and H (Figs. 1, 2, 3, 4, 5). 

When both GC and WL data are included, differ- 
ent combinations of the potentials $ and ^> are probed 
from independent measurements of p g S g and p m S m . It 
helps break the degeneracies (in particular with the 
bias parameters) significantly. When adding the cross- 
correlations of GC and WL, further information about 
the time evolution of the perturbations at multiple red- 
shift bins is provided. In fact, as we will see, the cross- 
correlation data puts the strongest constraints on the MG 
parameters and breaks the degeneracies best. 



A. f(R) models parametrized by Bo 

Bo is related to the Compton wavelenght of the scalar 
degree of freedom (DoF). Below this transition scale, the 
growth of structure is modified due to the presense of 
the fifth force. It should be kept in mind when one wants 
to explain the correlations of Bo with other parameters. 
Fig. 1 shows the normalized (left) and un-normalized 
(right) correlations between Bo (#7) and the rest of the 
model parameters. 

As expected, the constraints (shown by scattered plots) 
are tighter and the correlations (the bar plots) are de- 
creasing when new data included. With all data, except 
for the bias parameters, the correlations (black bars) are 
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FIG. 1. Correlation between -Bo (#7) and other parameters for different combinations of data, normalized (left) and un- 
normalized (right). 



about ~ 10% or less when compared to the parameter 
constraints (the \/Ub~b x Cjj values shown with black 
squares) . 

flbh 2 plays an important role at the recombination, so 
that the height and separation of the CMB peaks are 
very sensitive to it. Hence, il^h 2 is well constrained by 
the CMB data and is not highly correlated with B . The 
correlation is significant only when the WL data is con- 
sidered without the GC data. In this case the constraints 
on all parameters, including Hf,h , are relatively tight but 
there is degeneracy between Bq and ilbh 2 , ft c h 2 and Hq. 
Such degeneracies are broken when GC data is included 
(Figs. 1-2). The CDM energy density is important af- 
ter the matter-radiation equality epoch 2 . CMB data is 
mainly sensitive to Sl c h 2 through the geometry of the 
universe. The SNe la data is also sensitive to it through 
the value of the Hubble parameter entering the distance 
measurements. Varying Q c h 2 will change the epochs of 
matter-radiation and matter-dark energy equalities, lead- 
ing to a different history of growth of perturbations. The 
CDM density is also important for the secondary CMB 
effects, such as the ISW effect and is an extra source of 
correlations between the MG parameter, Bq, and il c h 2 . 
Hq enters in both the distance measurements and the 
evolution equations of matter perturbations, Eq. (8). As 
a result, the MG parameters are highly degenerate with 
H after including the large scale structure data. With 
the cross-correlation data, the background evolution is 
discriminated from the matter density perturbations and 
the degeneracy with Hq is significantly reduced. 

The power spectrum of the primordial perturbations is 



2 Q c is not directly measured and is inferred from the measuements 
of total matter or baryonic matter densities. 



related to n s and A s parameters as 

oc A s k n >~ A , (19) 

so that they determine the initial values of scalar metric 
perturbations, as well as the matter density perturba- 
tions. The impact of a modification to growth after re- 
combination can, in principle, be degenerate with varying 
n s and A s . The main impact of the optical depth, r (not 
to be confused with conformal time), is on the CMB, in- 
cluding the ISW term, which can be degenerate with evo- 
lution of the density perturbations. With the GC and WL 
data, this degeneracy still exists while the constraints on 
the parameters become tighter. We therefore see higher 
correlation in Fig. 1 (left), r is also correlated with n s 
so that tight constraints on n s can improve errors on r. 
Again, with the cross-correlation data, the time evolu- 
tion of perturbations is much better constrained and the 
degeneracy is broken. 

The EoS of DE, w, like other background parameters, 
affects the growth by changing the expansion history. The 
dark energy density, and consequently w, influences the 
recent expansion history and determines the ISW signal. 
However, since most of the information about w is coming 
from SNe la, the main role of the large scale structure 
data is to break degeneracies with the MG parameters 
(Fig. 2). 

As can be seen from Fig. 1 and later, the bias parame- 
ters would have the highest degeneracies with the MG pa- 
rameters. Observations of galaxies (baryonic matter) are 
related to the matter density perturbations via the bias 
parameter. GC data alone is neither able to constrain the 
bias parameters, nor it can break their degeneracies with 
Bq. Eq. (18) shows that increasing Bq and decreasing 
the bias have the same impact on the growth of perturba- 
tions. This negative correlation, however, would vanish at 
high redshifts due to less observable galaxies. When WL 
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FIG. 2. The 68% confidence level contours for the basic parameters and Bq. CMB+SNe data is included in all the cases. 



data is added, the sensitivity to the bias parameters, spe- 
cially at high redshifts, is increased. However, if we look 
at the uncorrelated plots, we see that the constraints on 
parameters are shrunk by almost two orders of magni- 
tude. Hence, the residual degeneracy, although relatively 
high, are not as important. The cross-correlation data is 
able to lower the uncertainties and the correlations by a 
factor of few. 

While the above discussion provides insight about how 
our model parameters could be correlated with the MG 
parameters, with the 1-parameter MG model we can 
mainly study the theories where there is an intrinsic tran- 
sition scale in the behavior of gravity due to a scalar de- 
gree of freedom with a fixed (f(R)) coupling to matter. In 
the following subsections, we investigate the degeneracies 
for models with more MG parameters. 



B. Sealer-tensor models parametrized with (c, s,n) 

Fig. 3 shows the correlation plots for c (top row), s 
(middle row) and n (bottom row). Here, c (#7) is a 
pre-factor controlling the overal amplitude of MG effects, 
while s (#8) and n (#9) control the redshift and scale 
of the MG effects, respectively. This parametrization is 
more flexible in capturing a MG feature, compared to 
the 1-parameter model in the previous section. However, 
the MG parameters themselves are highly correlated. The 
high correlations limit how well the MG parameters can 
be constrained, and one sees that adding the WL data 



does not improve the constraints significantly. 

From the left panels of Fig. 3, we see that c has a neg- 
ative correlation with both s and n. An increase of c will 
enhance the growth on all scales. Increasing n decreases 
the enhancement of growth on larger scales and pushes 
the MG effects to smaller scales. Hence, the impact of c 
on the linear growth is opposite to that of n. Increasing s 
pushes the MG effects to more recent times (smaller red- 
shifts). Data is able to discriminate between c and s from 
the measurement of the growth at multiple redshift bins 
but c and n are producing highly degenerate effects, s and 
n, on the other hand, are positively correlated because 
increasing s pushes the MG effects to smaller redshifts. 
For the WL data, as the dominant observable, this effect 
is degenerate with increasing n, which pushes the MG 
effects to smaller scales. This degeneracy is an intrinsic 
property of the WL potential, as described in [38]. 

Fig. 3 shows that degeneracies of the MG parame- 
ters with other parameters are broken significantly af- 
ter including the WL data. In particular, degeneracies 
with the bias parameters are reduced so that, Cjj — 
Cij/\J{Ca x Cjj) < 0.1 in all cases. One can see that 
there is a sweet spot in redshift where the bias parame- 
ters are best constrained. This spot in the middle redshift 
range is where the experiments are most sensitive to MG 
features [38]. 

As we see, adding new MG parameters provides more 
flexibility in detecting MG effects but the MG parameters 
will be highly correlated. One suspects that adding even 
more MG parameters to the model might just weaken 
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the constraints, without providing any new information. 
In the next section, we consider a model with more MG 
parameters and test how redundant the extra parameters 
would be. 



C. General model parametrized with BZ 
parameters 

Fig. 4 shows the correlations for log(Af) (first row), 
log(A|) (second row), fli (third row) and s (last row). 
Here, (/3i,/?2) are the analogues of c in the (c, s, n) model, 
while (Ai,A2) are similar to n, i.e. controlling the transi- 
tion scale where the MG effects become important, s in 
both models has the same role. 

The main reason for including this model is to show 
the redundancy of the extra MG parameters introduced. 
Comparing the plots in the Figs. 3 and 4 shows that 
there are more degeneracies in this model, specially with 
the bias parameters. We notice that Ai and A2 are con- 
strained very similarly by data. It has been shown that 
fi (and hence, Ai) should be better constrained than 7 
(hence, A 2 ) by the observables we use here [38]. But, we 
see in Fig. 4 that there is almost a perfect correlation 
between Ai and A 2 . This correlation limits how well data 
constrains Ai, via the uncertainties of A2 measurement. 
A's also have a high positive correlation with s (#11) 
due to the same reason that explained the correlation 
between n and s in the (c, s,n) model. Similarly, /3i and 
P2 are also highly correlated and we did not include the 
correlation plots for /3 2 . 

As before, larger s would lead to the enhancement of 
the growth starting at a larger redshift and accumulating 
over time. This mimics the effect of a stronger fifth force 
due to an increase of /3's. Since s is a common parameter 
to both the BZ model of this section and the (c, s,n) 
model of the previous section, with the same fiducial 
value, we can get an estimate of how much the constraints 
on the parameters are diluted and how much extra degen- 
eracy is introduced. Comparing the un-normalized plots 
for s in Figs. 3 and 4 confirms that both the errors and 
the correlations are higher by at least one order of mag- 
nitude for the BZ parametrization. 



D. Pixelated MG functions 

In Fig. 5, the correlation plots are shown for the 
middle-k (0.07 < k < 0.14 h/Mpc), low-z pixels of ^ 
(top row) and 7 (bottom row). We choose these pixels 
since they cover the scales and redshifts where LSST is 
most sensitive to the modified growth [38]. 

In the previous models, the MG parameters were in- 
troduced in a way that they would affect the growth at 
all scales and redshifts. This inevitably produces corre- 
lations between the MG parameters. On the other hand, 
with the MG pixels, the modified growth is confined to 
the ranges covered by each pixel so that correlations are 



intrinsically lower. In particular, the fi pixels are less cor- 
related with the 7 pixels compared with the correlations 
in the BZ model. As a result, the errors on the fi pixels 
(#8-13) are obviously smaller than (and are not limited 
by) those of the 7 pixels (#14-19). 

The other advantage of working with the MG pixels 
is that one can isolate some of the physical processes 
that affect the growth at particular scales or redshifts. 
For example, we see from Fig. 5 that there is a relatively 
smaller error on the first high-z 7 pixel (# 15) compared 
to the error on the other two pixels (# 17,13). This is 
due to the ISW effect which strongly constrains the 7 
function at large scales [38] which is covered by the first 
high-z 7 pixel. 

One can also look at the correlations between the MG 
pixels and w bins (Fig. 5). There is a high correlation 
between the MG pixels and the high-z w bin (#7), which 
is positive (negative) for /i (7). This is reasonable as in- 
creasing w would change the onset of the acceleration 
phase to a later time and hence, more time for the growth 
of perturbations. On the other hand, the low-z w bin (#6) 
has a negligible degeneracy with the MG pixels, again, in 
accordance with what we expect. Since most of the low-z 
information on w is coming from SNe la data and, the 
low-z w bins are not correlated with the MG pixels. 

The plots also show a general trend of negative correla- 
tions between the high-z pixels with their corresponding 
low-z pixels. It makes sense as any change in the high-z 
(low-z) pixel value should be counteracted with a change 
in the corresponding low-z (high-z) pixel values. As for 
MG parameters in the previous models, there are also 
high correlations between pixels and basic parameters, 
like H or A s (#20), and the bias parameters. 



IV. OTHER DEGENERATE EFFECTS 

In this subsection, we take the 3-parameter (c, s,n) 
model and generalize our analysis by including the WL 
systematics and considering neutrinos to be massive. We 
would like to know how these additional effects change 
the correlations between the MG and other parameters. 
We first include the WL systematics and, in the next 
step, the massive neutrinos. 

The systematics we consider here are the photo- z er- 
rors and some of the errors in the measurement of the 
point spread function (PSF). These errors are modeled 
in [27, 28] and we use their parametrization for our Fisher 
analysis (see [38] for details). 

Massive neutrinos are taken here as a part of the dark 
matter (DM) in the universe and their total mass is con- 
sidered as a parameter in the Fisher analysis, with a fidu- 
cial value of £m„ = 0.05 eV. 

Fig. 6 shows the correlations for c (top row), s (middle 
row) and n (bottom row) in four cases: without any ad- 
ditional effects, with systematics, with massive neutrinos 
and with both systematics and massive neutrinos. For 
these plots we used the full data combination. The to- 
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tal mass of neutrinos is considered as the last parameter 
(#21) in the plots. 

Systematics do not change the parameter constraints 
and their correlations considerably. With massive neu- 
trinos, however, the correlations of MG parameters with 
each other, with the bias parameters and some of the ba- 
sic parameters increases significantly. Massive neutrinos, 



as part of the DM, would freely stream and slow down the 
growth on the free-streaming scales. The heavier they are, 
the larger portion of DM they make, leaving less CDM 
in the universe. This slows down the growth and, hence, 
massive neutrinos are negatively correlated with c. Heav- 
ier neutrinos would have a shorter free-streaming scale. 
This is roughly the same as the effect of increasing s or 
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n, thus, explaining their positive correlation. 

V. SUMMARY 

In this paper, we studied the degeneracies between dif- 
ferent sets of MG parameters with the other parame- 
ters that are usually considered in models of large scale 
structure formation. We used a combination of the CMB, 
SNe la and large scale structure experiments, performed 
a Fisher analysis of the parameters and computed their 
correlation matrix. 

We considered three models with parametric forms 
of the MG functions. We saw in Section III that in 
the 3-parameter MG model, degeneracies were broken 
to a good extent but there were relatively high correla- 
tions between the MG parameters themselves. It was also 
demonstrated that increasing the number of the MG pa- 
rameters would introduce redundancies and weaken the 
constraints. Finally, with pixelization of the MG func- 
tions, we were able to study the degeneracies at different 
scales and redshifts. 

Among the model parameters, the bias parameters had 
the highest correlation with the MG parameters, due to 
their direct involvement in the equation for the growth of 
perturbations, Eq. (8). It seems that for the 3-parameter 



model, it was possible to break such degeneracies better 
without weakening the overall constraints on the param- 
eters. The background parameters were degenerate with 
MG mainly through the expansion history and its im- 
pact on the growth time. With the GC and WL data 
combined, one probes two different combinations of the 
scalar potentials. This leads to tight constraints and rela- 
tively low degeneracies. Among the observables we used, 
WL provides most of the information about the growth 
at low redshifts. The intrinsic degeneracy between the 
scale and the redshift in the WL kernel limits its abil- 
ity in probing the time evolution of perturbations. This 
leads to some residual degeneracies between the MG and 
the basic parameters. In such cases, the cross-correlation 
of WL and GC data at multiple tomographic redshift 
bins will significantly break the degeneracies. In fact, the 
cross-correlation power spectra are the most informative 
observables for probing the growth of structure. 

We have also considered other physical effects that 
could be degenerate with a modification to gravity. The 
WL systematics that we considered in our analysis did 
not change the correlations while the presence of massive 
neutrinos could have an important impact. It would be 
interesting to investigate how including new type of ob- 
servables, such as redshift space distortion (RSD) [8, 44], 
would resolve such degeneracies. 
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